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THE FIELD EQUATIONS 



The colliding wave problem has been described in terms of two approach- 
ing waves in regions II and III, in a background region I, which is here 
taken to be flat. According to the work of Penrose (1980), the initial data 
are well set, so that a unique solution exists in the interaction region IV 
at least in the neighbourhood of the boundaries of regions II and III. It is 
therefore necessary first to state the relevant field equations in the interac- 
tion region, and then to attempt to solve them subject to the appropriate 
boundary conditions. This formal approach is not always achievable in 
practice, as analytic solutions in region IV are not easily obtained. It is 
often more appropriate to start with the field equations for region IV and, 
once a solution is found, subsequently to consider the particular initial 
waves that would give rise to it. In either case, it is necessary to formulate 
the appropriate field equations for the interaction region. The purpose of 
this chapter is to derive such equations. 

Following Szekeres (1972) and Griffiths (1976&), the field equations 
will be derived here using Newman-Penrose techniques. 

6.1 The coordinate system 

The approaching waves have plane symmetry. On collision they will warp 
and scatter each other. In the interaction region they will no longer retain 
the full plane symmetry described by (4.15), though it may be expected 
that they will retain a two-parameter symmetry group of motions gener- 
ated by the Killing vectors 

6 = dx, 6 = dy (6.1) 

We expect the interaction region to retain plane symmetry in this sense. 
We therefore initially assume that the Killing vectors and ^2 exist 
throughout the space-time, and that the metric therefore has no explicit 
dependence on the coordinates x and y in the interaction region. Ul- 
timately of course, this assumption is only justified if it leads to exact 
solutions. 

At each point of the manifold there exist just two null directions 
orthogonal to the planes spanned by and dy. We may now align the 
two vectors l'^ and n'^ with these two directions. Since the congruences 
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in these directions are hypersurface orthogonal, the vectors and are 
proportional to gradients. It is therefore possible to put 

= A-^u,^, = B-^v,i, (6.2) 

where A and B are independent of x and y. 

We may now choose u and v as suitable null coordinates in the in- 
teraction region. At this point, these are arbitrary functions of the co- 
ordinates with the same labels in regions II and III continued into the 
interaction region. It is therefore appropriate to assume that u > and 
V > m region IV. Labelling the coordinates as (u, v, x, y) = (a;°, x^, 
x^, x^), the null tetrad may now be chosen to be 

l^ = A-^5l, l" = {0,B,Y^,Y^) 

B-'Sl, n^^{A,0,X^X') (6.3) 



where, putting i = 2,3, 

X' = X'{u,v), Y'^Y\u,v), C = C{u,v). (6.4) 

The tetrad and coordinates are now defined up to the following trans- 
formations: 

(1 ) Scale transformations or boosts: 

l^" = (j)l^', n^" = (t)-^ni' , (j) = (j){u,v), (6.5) 
under which the scale functions A and B transform as 

A' = (i)-^A, B' = 

(2) Spatial rotations or spins: 

m^" = e'^m^, C = C{u, v). (6.6) 
This is a rotation in a space-like 2-surface and transforms the ^* by 

e = e'^'e- 

(3) Relabelling of null hypersurfaces: 

u' = f{u), v' = g{v), (6.7) 
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which induces the transformation in the scale functions 



A' = l^A , B' = g,,B. 

(4) Spatial coordinate transformations: 

u' = u, v' = v, x'' = x' + h\u,v), (6.8) 
which transforms the X* and components by 

X'' = X' + h%A , Y'' = Y' + h%B. 

(5) Linear coordinate transformations: 

a;*' = a''jX\ a*^- = consts. (6.9) 

6.2 The derivation of the field equations 

With the tetrad defined as above, it may be noticed that all components 
are at most functions of u and v only. It follows that the differential 
operators, when applied to spin coeflBcients, take the form 

D = Bi;, A = Ai^, 6 = 0. (6.10) 

The incoming waves have no dependence on the coordinates x and y, and 
so no such dependence is expected in the interaction region. 

The metric equations, which are the commutation relations (2.15) 
applied to the coordinates, immediately take the form 

K = 1/ = 0, P — P, P = P 
f=7r = 2p = 2a 
DA= -{e + e)A 

AS = (7 + 7)5 (6.11) 
AY' - DX' = (7 + 7)y* + (e + e)X' - 4aC - 4af 
De = {p + e- e)e + aC 
Af = -(// + 7-7)f -Af. 



The first of these equations simply expresses the condition, consistent 
with (6.2), that the congruences tangent to l^^ and n^^ are geodesic and 
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twist- free. The third and fourth metric equations (6.11) enable the real 
parts of e and 7 to be defined by 

e + e = -S(logA),„ , 7 + 7 = A(logS),,. (6.12) 

It can thus be seen that the two null congruences do not necessarily have 
affine paramctrization. 

Before proceeding to the field equations, it is convenient to introduce 
modified 'scale-invariant' spin coefficients 
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It is also convenient to introduce the scale-invariant function M, defined 

by 

M = log vis. (6.14) 

The metric equations (6.11) and the Ricci identities (2.14) now only 
involve quantities that are invariant with respect to the scale transforma- 
tions (6.5). Together they give the following: 

y°% - = -4(a°f + a°e)e-^/2 (6.15a) 
e% = (p°+ii?°)e + aT (6.156) 
e,u = -{l^° -^G°)C -\e (6.15c) 

p° ^ = -2p°/i° - Aa°a° - - 3A° (6.15e) 
= 2p>° + 4a°a° + + 3A° (6.15/) 

= - A°A° - ^°M,„ - (6.15^) 

(7% = a°(2p° - M° + 2iE°) + % (6.15/1) 

(7% = -c7°(//° - 2iG°) - p°A° - 4a°2 - $°2 (6.15z) 

A% = A°(p° - 2zE°) + p,°a° + 4a°2 - (6.15j) 

A% = -A°(2/i° + + 2zG°) - -^l (6.15/c) 
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„ = Q;°(3p° - |M,„ - iE°) + a°a° + $°o 
-a°{Sii° + |M,^ + iG°) - a°X° - 
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(6.150 
(6.15m) 
(6.15n) 
(6.15o) 
(6.15p) 
(6.15?) 



Following the method of Szekeres, it is convenient at this point to 

put 

^2 = -i=e(^-^)/2(cosh|W^ + zsinh|W^) 

^3 ^ -i.e(^+^)/2(sinhiW^ + zcoshiW^) ^^'^^^ 

where U, V and W are functions of u and v. Using equations (6.156, c) 
some of the scale-invariant spin coefficients can then be expressed explic- 
itly as 



E° = -lVySinhW G° 
a° = \iWy -\Vy cosh W X° 



-|KsinhW^ (6.17) 



where the suffices u and v refer to partial derivatives taken with respect 
to these variables. From now on, the commas on metric functions will be 
omitted and derivatives will be denoted only by suffices. 

The method of procedure from this point depends heavily on the 
form of the Ricci tensor. For colliding gravitational and electromagnetic 
waves it will be shown that ^»oi and $21 are both zero. In this case it 
follows from (6.15/,m) that, since q;° = in the initial regions II and III, 
it is necessary that a° = everywhere. 

By contrast, for colliding neutrino fields (Griffiths 1976a) which will 
only be considered in Section 20.5, the components $01 and $21 are es- 
sentially non-zero and the method of integrating the field equations differs 
significantly. In this, and similar cases, it is clearly not possible to put 
a° = everywhere. Apart from section 20.5, however, the methods for 
integrating the field equations in these cases will not be pursued. 

In this book the emphasis is on gravitational and electromagnetic 
waves, so only the Einstein-Maxwell equations will be derived explicitly in 
this section. Einstein's vacuum equations are obtained simply by putting 
the electromagnetic field components equal to zero. 
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6.3 The Einstein and Einstein— MgLxwell equations 

In order to consider the collision of electromagnetic waves, it is convenient 
first to define the scale-invariant quantities 

$o = $o5"\ $1 = $^ = $2^"^ (6.18) 

With the conditions (6.11), Maxwell's equations (2.18) now become 

It is immediately clear from the first and last of these equations 
that, since $i is zero in regions II and III, it must also be zero in the 
interaction region. It follows that ^»oi $21 both zero for colliding 
electromagnetic waves, as they obviously are for pure gravitational waves. 
With this condition, equations (6.15 /, m) imply that q;° is also zero in 
region IV, since it is similarly zero in regions II and III. 

Equation (6.15a) now implies that, since a = 0, it is possible to use 
a spatial coordinate transformation (6.8) to make both X* and zero 
simultaneously. With this transformation, the metric now takes the form 

ds^ = 2e-^dudv 

TT O ,^0 (6-20) 

- e-^(e^ cosh Wdx^ - 2 sinh Wdxdy + e'^ cosh Wdy^). 

This will be referred to as the Szekeres line element. 

The only non-zero spin coefficients are now p, cr, e, A and 7. 
Apart from the freedom associated with the scale transformation (6.5), 
these are now expressed explicitly in terms of the metric functions by 
equations (6.17) and (6.12). 

Using (6.17) the remaining components of Maxwell's equations (6.19) 
now become 

^°2,v = Wy+ iVy sinh W)^l - \ {iWu + Vu cosh W)% ^ ^ 

1 1 (6.21) 

^o,u =\{Uu- iVuSmhW)^l + \{iWy - VyCoshW)^l. 

The expressions (6.17) may also be substituted back into equations 
(6.15 d,ej,g,ij,n) to give 

Uuv = UuUy (6.22a) 
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2U^^ = + W"„2 + y„ 2 cosh^ W - 2UyMy + 4$°§° (6.226) 
2Uuu = + Wj" + V^^ cosh^ W - 2UuMu + (6.22c) 



+ 2{^°o^ + ^^°Q)sechW (6.22d) 
2Wuv = UyWy + UyWu + 2VuVy suih W cosh W 

+ 2i{^l^l-^m) (6.22e) 

2Muv = -UuUy + WuWy + VuVy cosh' W. (6.22/) 



These are Einstein's field equations written in component form. 

The remaining equations (6.15 h,o,p,q, k) give the scale-invariant 
components of the Weyl tensor in the form 

^rg = -|((Kz; - UyV^ + M^Vy) cosh + 2VyWy siiihW) 

+ li{Wyy - UyWy + MyWy " ' COSh W Sluh W) 

= 

% = ^Muy - li{VuWy - VyWu) cosh W (6.23) 

n = -l {{Vuu - UuVu + MJ/u) cosh W + 2VuWu sinh W) 
- \i{Wuu - UuWu + MuWu - cosh sinh W^). 

It follows from this that it is only the Coulomb component \E'2 that ap- 
pears as an interaction term when the gravitational waves ^'o ^ind 
interact. 

6.4 Integrating the field equations 

Having derived field equations for region IV of the colliding plane wave 
situation, it is initially necessary to question the possible existence of solu- 
tions beyond the initial boundaries with regions II and III. This question 
has been considered by Yurtsever (1989a). In an appendix to this paper, 
he has given a formal proof of the global existence and uniqueness of so- 
lutions in the interaction region as far as the focusing singularity, subject 
to the boundary conditions which will be discussed in the next chapter. 

With this question settled, it is appropriate here to make a few simple 
observations about the basic approach to the integration of the above field 
equations. 

To start with, it can be seen that equation (6.22a) can immediately 
be integrated to give 

= f{u)+g{v) (6.24) 
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where f{u) and g{v) are arbitrary functions that, at this stage, may be 
assumed to be piecewise C^. 

It may also initially be observed that the Szekeres line element (6.20) 
is similar to that for the plane wave metric in Rosen form (4.13) that is 
appropriate for regions II and III, except that the metric functions are 
here functions of both u and v. In fact it is possible to use (6.20) to 
describe the space-time in all four regions, where U, V, W and M are 
functions of u and v in region IV, are functions of u only in region II, are 
functions of v only in region III, and are constants that can be transformed 
to be zero in region I. Regions II and III contain approaching plane waves 
that only differ from those described in Chapter 4 by a simple coordinate 
transformation (6.7) in u or v. 

Using this approach, it is appropriate to choose / and g such that 

/ = i for « < ^ ^ 

. (6.25) 
g = ^ for V <0. 

The approaching waves are thus partially described by the functions 
f{u) in region II, and g{v) in region III. Of the remaining functions, M 
can be removed in regions II and III by the transformation (6.7), and V 
and W must satisfy either (6.22c) or (6.226). The approaching waves are 
thus characterized by two independent functions in each region, which 
together describe the amplitude and polarization of the two approaching 
waves. 

Returning to consider the full field equations in region IV, it may 
be observed that equations (6.22(i, e), together with Maxwell's equations 
(6.21), are the integrability conditions for the field equations (6.22 6, c, /). 
Thus if functions U, V, W, $2 ^^'^ ^0 found to satisfy (6.22a, d, e) 
and (6.21), then a function M that satisfies the remaining equations 
(6.226, c, /) automatically exists. It should not be assumed that this func- 
tion should be zero. It is not easy to find analytic solutions of equations 
(6.22(i, e) and (6.21) and, once a solution is found, a particular M{u,v) 
can be obtained by integrating (6.226, c). In view of the required conti- 
nuity across the boundaries fx = and v = 0, non-zero functions M{u) 
and M{v) will usually occur in regions II and III respectively. Normally, 
it will be possible to reduce these to zero only in these regions, and then 
only for purposes of interpreting the solution obtained. 

It may be seen from this discussion that, in order to integrate the 
field equations, attention is focused on equations (6.22 d, e) and (6.21) 
or, for colliding gravitational waves, just on (6.22 d, e). The emphasis is 
therefore on these equations as we seek to derive the exact solutions given 
in subsequent chapters. 



